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Attractive interactions between fermions can 
produce a superfluid ground state, in which pairs 
of up and down spins swirl together in a coordi- 
nated, coherent dance. How is this dance affected 
by an imbalance in the population of up and down 
fermions? Do the extra fermions stand on the 
sides, or do they disrupt the dance? The most 
intriguing possibility is the formation of a modu- 
lated superfluid state, known as an LO phase, in 
which the excess fermions self-organize into do- 
main walls where the pairing amplitude changes 
sign. Despite fifty years of theoretical and exper- 
imental work, there has so far been no direct ob- 
servation of an LO phase. Here we propose an ex- 
periment in which two fermion clouds, prepared 
with unequal population imbalances, are allowed 
to expand and interfere. A zipper pattern in the 
interference fringes is unequivocal evidence of LO 
physics. Furthermore, because the experiment 
is resolved in time and in two spatial directions, 
we expect an observable signature even at finite 
temperatures (when thermal fluctuations destroy 
long-range LO order averaged over time). 

The study of modulated superfluidity has its roots in 
the context of a superconductor in a parallel magnetic 
field. The inter-electron attraction favors a superfluid 
(SF) state consisting of pairs of up- and down-spin elec- 
trons, whereas the field favors a polarized Fermi liq- 
uid (FL) state with a lower Zeeman energy. The sim- 
plest theories predict a first-order SF-FL transition^. 
However, the competition between pairing and polariza- 
tion can produce far more subtle physics - an interme- 
diate Larkin-Ovchinnikov (LO) state*™!, as depicted in 
Fig. 1. As the field is increased beyond h c \ it forces ex- 
cess fermions into the superfluid in the form of domain 
walls - topological defects - at which the order param- 
eter changes sign between positive and negative values. 
(This is analogous to how a perpendicular field forces 
vortices into a superconducting film.) The wavelength of 
these modulations decreases with increasing field and ul- 
timately the system gives way to a polarized Fermi liquid, 
shown in Fig. [I] 

As a fascinating example of self-organized quantum 
matter, LO states have long been sought after in su- 
perconductors, in ultracold atomic gases^, and even in 
neutron stars^. In the context of cold atoms, there have 
been proposals to detect LO states through (a) modula- 
tions of the polarization in real space; (b) peaks in the 
pair momentum distribution at the modulation wavevec- 
tor; (c) shadow features in the single-particle momentum 
distribution; and (d) Andreev bound states in the density 



of states. A recent experiment on an array of tubes found 
density profiles in agreement with Bethe ansatz calcula- 
tions, 9 which predict power-law LO correlations at zero 
temperature. However, direct evidence of the modula- 
tions of the order parameter - the defining feature of an 
LO state - is still lacking. 
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FIG. 1: (Top) Depictions of a fully paired superfluid, an LO 
state with excess fermions in domain walls, and a polarized 
Fermi liquid. (Center) Pairing amplitude A and magnetiza- 
tion m as a function of the Zeeman field h, which is the differ- 
ence between the chemical potentials of up and down spins. 
The LO phase exists in a field range h c \ < h < h C 2- (Bottom) 
Real space behavior of A(x) and m(x) in each phase. 

Experimental Proposal: In this Letter, we propose an 
interference experiment in which two isolated superflu- 
ids expand into each other, as illustrated in Fig. [2] The 
ideal situation is shown in the lower panel of this fig- 
ure, where one layer is a uniform fully-paired superfluid, 
which serves as a reference phase, and the other layer is 
a modulated (LO) superfluid. The resulting interference 
pattern is directly sensitive to real-space modulations of 
the order parameter, and should provide an unequivocal 
signature of the elusive LO phase. 

LO st ates ha ve been predicted to exist in various 
situation d 8 * 11 * 12 *. The most likely of these to be realized 
in the near future is an array of tubes with a small inter- 
tube coupling EHUD This quasi-one-dimensional geometry 
provides good Fermi surface nesting at the LO wavevec- 
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FIG. 2: Proposed experiment and interferometric LO signa- 
ture. An optical lattice is used to prepare two independent 
layers. These clouds are subsequently allowed to expand and 
interfere and are then imaged along one of the in-plane direc- 
tions (here y). If both clouds are in the fully paired SF phase, 
the interference pattern is similar to the familiar double-slit 
experiment (top). However, if one of the layers contains an 
LO state (the change in sign of A(x) is indicated by color), 
the interference pattern contains sharp dislocations of the in- 
terference fringes (bottom). 



tor #lo = fcpt — kpi, together with Josephson coupling 
between the order parameter in adjacent tubes which is 
necessary to stabilize true long-range order. 

Thus, we propose the following experiment. Two 
species of fermions (referred to as t and |) are loaded into 
an optical trap. The cloud is separated into two indepen- 
dent quasi-2D "pancake" layers using an optical lattice 
with a wide spacing in the vertical direction z (created 
with two laser beams intersecting at a shallow angle) 
The two layers are caused to have different population 
imbalances, such that n* op ^ n* op but n^ ottom = n J ottom . 
This may happen by chance due to natural number fluc- 
tuations; alternatively, one can induce hyperfine tran- 
sitions at the beat frequency between two laser beams, 
using a shallow-angle interference technique with appro- 
priate phase shifts to address each layer separately. An 
additional optical lattice is further used to create a 2D 
array of weakly-coupled tubes, conducive to the forma- 
tion of an LO state. The resulting geometry is shown in 
Fig-i 
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FIG. 3: The fermion gas is confined in a harmonic trap. An 
optical lattice with a large spacing in the ^-direction is used 
to separate the gas into two independent quasi- 2D layers. A 
second optical lattice in the ^/-direction cuts each layer into a 
series of weakly coupled tubes - the optimal geometry for LO 
physics. The trap and lattices are turned off abruptly, allow- 
ing the two layers to expand and interfere with one another. 
Finally, the shadow of a probe laser in the ^/-direction gives 
the interference pattern projected onto the x-z plane. 

Once the gases are allowed to equilibrate for sufficient 
time, the interactions in the system are quickly ramped 
to the BEC side of the Feshbach resonance 17 to "freeze" 
or "project" the pair wavefunction into a boson wave- 
function, so that from this point onwards the pairs move 
as independent bosons (instead of disintegrating into 
fermions). Then, the confining potentials are abruptly 
turned off. As the clouds expand into one another, they 
interfere and form a 3D matter wave interference pattern. 
The projection of this interference pattern onto the x-z 
plane can be measured by absorption imaging along the 
y direction. Any LO phase modulation features will be 
captured in these projected interference patterns. 

Interference between coupled tubes: We now discuss an- 
alytically the salient features of the interference patterns 
of a 2D array of coupled tubes. We begin by consider- 
ing two layers, each containing N coupled tubes, with 
separation d in the z-direction. After ramping up the in- 
teraction to produce a molecular BEC (of fermion pairs) , 
the wavefunction is 

A(x,y, Z ) = e -(*-dmV^ST A T {x)e -(y-anf/2*l 

n=l 

top layer 

N 

+ e -(*+d/2) a /2«r; J2 A* (x)e-( y - an ^' 2 °Xl) 

n=l 

v v ' 

bottom layer 

where a is the separation between in-plane tubes and 
<jy and cr z are the Gaussian confinements in the respec- 
tive directions. The wavefunction in the n th tube is de- 
noted by A^(x) in the top layer and by A^(x) in the 
bottom layer. When the trap and lattices are switched 
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FIG. 4: Interference patterns for three different configura- 
tions: fully paired SF state (top), locked LO state (middle), 
and unlocked LO state (bottom). In each case we consider 
a 2 x 5 array of coupled tubes as shown in Fig. [3] using a 
fully paired SF as a reference phase in the bottom layer. The 
pairing amplitude A(x) in each of the five top layer tubes is 
shown to the left of the interference patternt. The locked LO 
states were taken to be A(x) = e~ x ^ 2cr sn(x/L\k), for k < 1 
with period L, where sn(x\k) is the Jacobi sine function. In 
the unlocked case we added a random displacement of the do- 
main walls A(x) = e~ x2/2(r2 sn((x+5)/L\k) where S e [-L,L]. 
This interference pattern still contains signatures of the LO 
phase even when the domain wall locations fluctuate between 
tubes. 



off, the clouds expand predominantly in the tightly con- 
fined directions (y and z). After a suitable time of flight 



t, the wavefunction is effectively Fourier- transformed in 
the y and z directions. That is, the final wavefunction 
A(x,y, z,t) is approximately proportional to the initial 
momentum distribution in the y and z directions, i.e., 
A(x,y,z,i) ~ A(x,k y ,k z ,t = 0) where y = tk y /m and 
z = tk z /m: 



A(x, k y , k z ) 



(2) 



X 



N 



ik z d/2 



n=l 
N 

n=l 

The 3D density of the cloud, after expansion, is given 
by I(x : k y: k z ) ~ |A(x, k y: k z )\ 2 . The imaging process 
measures the integrated density along the y direction, 
I(x,k z ) ~ / dky\A(x,k y ,k z )\ 2 : 

N N 

I(x,k z ) ~ e~^ 2 J2 E e-("- m ) 2 ° 2 / 4 ^ (3) 

n—1 m=l 

x [AZ(x)Al(x) + AB( x ) A B(x) 
+ 2Al(x)A^(x) cosk z d] . 

We can consider the behavior of the above interfer- 
ence formula in its two limits: widely separated tubes 
(a/a y — » oo ) and overlapping tubes (a/a y —> 0). In the 
overlapping limit, the total interference pattern reduces 
to that of two isolated 2D layers. In the limit of widely 
separated tubes, similar to the proposed experiment, the 
total interference pattern reduces to the sum of the inter- 
ference patterns from adjacent tubes in the two layers. 

The above analysis involves just two layers, which are 
necessary and sufficient for generating interference. In- 
troducing more layers would complicate the analysis, and 
reduce the visibility of the interference pattern; neverthe- 
less, there would still be observable effects of the order 

Of 1/v^layers- 

Figure [4] illustrates the projected interference pattern, 
described by Eq. |4j for three configurations of the up- 
per layer: a uniform SF; an LO phase with domain walls 
locked between tubes; and an LO phase with domain 
walls fluctuating between tubes. In each case we assume 
that the lower layer has been prepared in a uniform SF 
state. We represent the LO pairing amplitude modula- 
tions using Jacobi sine functions sn(x\k]^ multiplied by 
Gaussian envelopes in the x, y, and z directions. The 
"zipper" pattern in the lower two panels is a clear sig- 
nature of oscillations of the relative phase between the 
SF and LO layers, in contrast to the straight interference 
fringes in the top panel. 

We have also obtained similar results from Bogoliubov- 
de Gennes (BdG) simulations. The BdG simulations con- 
verge to "locked LO" states even when the intertube cou- 
pling is small. Thermal phase fluctuations beyond BdG 
may be expected to produce domain wall fluctuations be- 
tween tubes, as in the lowest panel of Fig. [4j but even 
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FIG. 5: Illustration of imaginary time worldlines of LO do- 
main walls and corresponding magnetizations with (bottom) 
and without (top) quantum fluctuations. Pluses and minuses 
correspond to the sign of the pairing amplitude A(x) in each 
region. In the bottom panel, the quantum fluctuations en- 
force a spatial profile for the magnetization that is less sharp, 
but which still clearly maintains the features of domain walls. 



then a zipper-like pattern is still visible provided that the 
fluctuations are not too severe. 

Discussion: In an experiment previously proposed by 
the cold atom group at Rice, an LO state is prepared in 
a 2D array of tubes (with a spacing of about 532nm), 
the interaction is ramped to the BEC side, the trap is 
switched off, and the final density profile is a measure 
of the initial longitudinal pair momentum distribution 
P(q x ) (along the tube axis). The function P(q x ) is ex- 
pected to have cusps at i^LO, but these features tend to 
get washed out by temperature and spatial inhomogene- 
ity, and their interpretation requires detailed quantitative 
analysis. The key advantage of our proposal is that the 
clouds expand in the transverse (z) direction, giving an 



interference pattern resolved along the x-axis, which pro- 
vides a much more direct way to probe the LO pairing 
amplitude modulations. 

The penetration of domain walls in an LO state re- 
sembles the penetration of vortices into a rotating BEC, 
and indeed, the zipper interference pattern we predict is 
similar to patterns that have been seen in experiments 
on vortices in 2D Bose systems^. Those experiments 
involved two pancake BECs separated by 3 /im. 

A remarkable feature of cold atom experiments is that 
control parameters (interaction, lattice depth, and trap 
depth) can be turned off very quickly, much faster than 
the typical timescale of domain wall movement. This 
allows us to take a snapshot of the wavefunction (resolved 
in real time), in a way not possible in condensed matter 
experiments. Thus, even above the critical temperature 
where thermal fluctuations destroy long-range order, it 
may still be possible to detect LO physics in the form of 
"temporary" domain walls! 

Quantum fluctuations of an LO state are more subtle. 
Their effect can be thought of as diffusion of the domain 
walls in imaginary time r; measurements are necessarily 
averaged over r. For isolated tubes, quantum fluctua- 
tions of the domain walls prevent long-range LO order, 
and there is only quasi-long-range order at zero tempera- 
ture; hence, the interference pattern will be washed out. 
This is why we recommend using sufficient coupling be- 
tween the tubes to stabilize long-range order, so that the 
pairing amplitude modulations remain even after averag- 
ing over quantum fluctuations (see Fig. [5|. 

In a larger perspective, interferometric techniques have 
proven to be powerful methods to detect the relative 
phase of the pair wavefunction in cupratesp^ruthenates, 
pnictides, and other unconventional superconductors.^ 
The interferometric signature we propose is different in 
that it involves a phase change in the center-of-mass mod- 
ulations of the order parameter that occur across the do- 
main walls and not in the relative phase of the pair. 
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